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Abstract 

We examine stochastic temperature fluctuations of the cosmic background ra- 
diation (CBR) arising via the Sachs- Wolfe effect from gravitational wave per- 
turbations produced in the early universe. These temperature fluctuations are 
described by an angular correlation function C{p(). A new (more concise and 
general) derivation of Ci^y) is given, and evaluated for inflationary-universe 
cosmologies. This yields standard results for angles 7 greater than a few de- 
grees, but new results for smaller angles, because we do not make standard 
long- wavelength approximations to the gravitational wave mode functions. 
The function C^y) may be expanded in a series of Legendre polynomials; we 
use numerical methods to compare the coefficients of the resulting expansion 
in our exact calculation with standard (approximate) results. We also report 
some progress towards finding a closed form expression for C(j). 

BACS numbers: 98.80.Cq, 98.80.C, 98.80.Es 
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I. INTRODUCTION 



Penzias and Wilson [0 discovered the Cosmic Background Radiation (CBR) in 1965. 
Since then researchers have studied the CBR using ground, balloon, rocket, and satellite 
based experiments 0||. The evidence indicates that this radiation is a remnant of an 
early hot phase of the universe, emitted when ionized hydrogen and electrons combined at 
a temperature of about 4000 K ||. In the simplest models this combination occurs at a 
redshift Z rs 1300, although it is also possible that the hydrogen was re- ionized as recently 
as redshift Z w 100 ||. In effect, the CBR is a picture of our universe when it was much 
smaller and hotter than it is today. 

The CBR has a thermal (blackbody) spectrum, and is remarkably isotropic and uniform. 
Only recently have experiments reliably detected perturbations away from perfect isotropy. 
Such perturbations are expected; in 1967 Sachs and Wolfe |5| showed how variations in 
the density of the cosmological fluid and gravitational wave perturbations result in CBR 
temperature fluctuations, even if the surface of last scattering was perfectly uniform in 
temperature. 

During the past several years, the Cosmic Background Explorer (COBE) satellite team 
has reported detailed measurements of the statistical properties of these temperature per- 
turbations |§. Analyzing the COBE data is subtle; it requires subtraction of the dipole 
and quadrupole moments arising from the Doppler shift due to the earth's peculiar velocity 
w.r.t. the cosmological fluid, and also the subtraction of infra-red and microwave emission 
from stars, dust clouds and gas within our own galaxy. In this paper, we assume that these 
contaminants have been removed from the data, and discuss only the perturbations of the 
CBR which are cosmological in origin. 

Additional measurements by other experimental groups |7|-[i4| have also reported per- 



turbations of the CBR over variety of angular scales. The range of angular scales covered 



by these different experiments is nicely illustrated in Figure lb of reference [15j ; the angular 
scales range from full sky coverage (180 degrees) down to angular scales less than 1/10 of 
a degree. These experiments are ongoing, and additional data should appear from these 
research groups over the next few years. 

For our purpose, the most useful statistical quantity determined by COBE, and the 
device frequently used to state and compare the results of the other experiments, is the 
sky-averaged angular correlation function 

C{u\v b ) = C(j) = (^K)^>)) sky . (1.1) 

In this formula, u a and v b are two unit-length spatial vectors, pointing out from the observer's 
location to points on the celestial sphere. The CBR temperature fluctuation in the direction 
u a away from the mean value of T is denoted 5T(u a ). The angle brackets, as used in 
individual experiments, refer to a uniform "sky average" over all points on the celestial 
sphere separated by angle 7, where cos 7 = u a v a . For this reason, the correlation function 
depends only on the angle 7, and not on the absolute position of the vectors u a and v b . It 
is convenient to expand this function in terms of Legendre polynomials: 

00 2/ -I- 1 

C(7) = £-^a^(cos7). (1.2) 
1=2 47r 
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The coefficients af are referred to as the multipole moments of the expansion. Note that the 
monopole term (I = 0) is absent; the dipole term (7 = 1) is generally removed from the data 
because it depends mostly upon the observer's peculiar velocity. The quantity measured by 
a given experiment is affected by the filtering properties of the optics and receivers, which 
determine the angular range over which the experiment is sensitive; the effect of this filtering 
is incorporated into the "weight function" Wi, which differs from experiment to experiment. 
These weight functions are shown in Figure lb of reference (Tj| for a number of different 
experiments; for the purposes of this paper we will consider an "ideal" experiment that is 
equally sensitive at all angular scales and has a weight function Wi — 1. 

With cosmological models concrete enough to make definite theoretical predictions, one 
may calculate the expected value of this correlation function. If the cosmological model is 
isotropic then the correlation function depends only on the angle 7 between the pair of 
observation points even before the "sky averaging" in (|1 . 1| ) is done; thus averaging is not 
necessary. It is important to note however that in the experimental case, the multipole 
moments af associated with the observed sky-averaged correlation function have definite 
measurable values, but given a specific theoretical model, these actual values are impossible 
to predict; they depend upon our location in universe, and additonally reflect the fact that 
our universe is a single realization of the statistical ensemble whose expected values may be 
determined theoretically. Hence, the quantity determined in this paper is an expectation 
value; we denote the associated expected multipole moments by (af). This ensemble average 
or expectation value is equal to a the uniform average over observers located at all spatial 
locations, with all possible choices of direction on the celestial sphere. Thus, in principle, 
one could directly compare the observed af with the expected values (af) by averaging 
observational data taken from regions of the universe that are currently not in causal contact, 
but this process would take many times the age of the universe to complete. Hence there 
remains a practical problem, that of constraining the cosmological models by comparing the 
expected multipole moments (af) with the observed multipoles af. This requires statistical 
analysis of the expected variance in af; this problem of cosmic variance will not be addressed 
here. One does expect however that since the number of independent degrees of freedom in 
the I'th multipole moment is 21 + 1, for a good cosmological model the observed af and the 
expected (af) should be very closely equal for large I. 

This paper considers the angular correlation function for models of the universe that 
pass through an early inflationary stage. More precisely, we consider models where the 
cosmological length scale (scale factor) undergoes a long period of exponential expansion, 
characterized by a constant, positive energy density and a constant, negative isotropic pres- 
sure of equal magnitude. Such cosmological models are attractive because they solve the 
horizon and flatness problems in a "natural" way [[L^,[ij]]. For this reason an enormous 
variety of mechanisms for inflation have been proposed during the past decade. 

Since the proposed inflationary models differ in significant ways, they make certain pre- 
dictions that are quantitatively very different. As an example, the perturbations in the 
CBR temperature that result from fluctuations in the matter-density are model-dependent 
because the matter content of these models is limited only by the imagination of the model- 
builder. Fortunately there are certain predictions of inflationary models that are independent 
of the details of the model. One significant example is the subject of this paper; the pertur- 
bations of the CBR temperature that result from the gravitational wave fluctuations. These 
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perturbations depend only on a single parameter: the energy-density during the period of 
exponential expansion. 

In the simplest perturbed FRW models, one may classify the perturbations which produce 
fluctuations of the CBR temperature as scalar, vector or tensor in nature. The complete 
angular correlation function C(j) is the sum of terms arising from each of these; one generally 
assumes that these add incoherently (or in quadrature). For a large class of "slow rollover" 
inflationary models, the expectation value of the angular correlation function resulting from 
scalar perturbations is 

(Note that the dipole moment has been removed.) This corresponds to an expected spectrum 
of coefficients (af) given by (af) = ^° 2 ^ . Only the overall amplitude of the correlation 
function, here determined by the expected value of the quadrupole moment (af), varies 
from model to model. The correlation function due to vector perturbations is typically 
very small and is neglected. In this paper we only consider the contribution to the angular 
correlation function from the tensor (or gravitational wave) perturbations; as we will shortly 
explain, these are entirely determined by the energy-density during the inflationary phase 
and are otherwise model- independent. 

There is a substantial body of research on this topic. In the following brief review we 
do not include much of the important work on the effects of scalar density fluctuations on 
the CBR, but principally discuss the work on CBR fluctuations induced by gravitational 
wave perturbations. The original discovery that cosmological expansion could create par- 
ticles is due to Parker |18| and Zel'dovich |19 |. However, they apparently assumed that 



the linearized gravitational wave equation would be conformally invariant and hence that 



no gravitons could be created. This oversight was corrected by Grishchuk |20| who showed 
that due to the lack of conformal invariance a period of rapid cosmological expansion could 
result in the non-adiabatic amplification of weak classical gravitational waves. The corre- 
sponding classical process for black holes (super-radiant scattering) implies the quantum 
effect (Hawking radiation). In similar fashion, Ford and Parker |TI| showed how one could 
systematically quantize the linearized gravitational field on an FRW background, and cal- 
culated the spectrum of gravitons created by the cosmological expansion. 

Starobinsky investigated this process in detail for inflationary cosmologies (but before 
the term "inflation" had been coined and before the advantages of such a period of 



expansion had been fully appreciated and explained [|I7]]) and found the power spectrum of 
gravitational radiation that would be left behind ||22j| . The quadrupole (/ = 2) and octupole 



(I = 3) anisotropies in the CBR induced by the resulting gravitational perturbations were 
later calculated by Rubakov, Sazhin and Veryaskin p3| . While the methods used and the 



interpretation of the results are entirely correct, this work suffers from technical errors. 
In particular, the octupole moment is correct but the quadrupole moment has the wrong 
value: the right-hand side of their equation (7) reads "2.4ey / 'M*", but the correct result is 
(AT /T )q Uadl . = 1.54ey/Mp 1 . Soon afterwards, similar results were published by Fabbri and 
Pollock pl| , who gave the first general formula for the I'th multipole moment in inflationary 
models. This work has a minor typographical error (the right-hand side of their equation (14) 
should be doubled) but otherwise their results are correct. About a year later, this work 
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was repeated and generalized for power-law inflation by Abbott and Wise f25fl , who also 
give a (now standard) correct formula for the I'th multipole moment. Shortly thereafter, 
Starobinsky EB] also published the results of an independent analysis, giving the same 
formula for the I'th multipole and correcting the errors in and ||24j| . This early work all 
considered the spatially flat k = FRW models; it was generalized to the k = ±1 cases by 
Abbott and Shaefer [^7], who systematically considered the CBR fluctuations induced by 
all three (scalar, vector, and tensor) types of perturbations to the k = 0, ±1 FRW metrics in 
inflationary models. The energy-density of the classical gravitational waves resulting from 
inflation was re-examined by Abbott and Harari who stressed the quantum-mechanical 
origin of this radiation, and by Allen |29[ who elucidated the first complete formula for 



the power-spectrum in gravitational radiation, and it's connection to the low- frequency 
instability (and peculiar infra-red behavior) of de Sitter space. As one consequence, Allen 
showed that the energy-density in gravitational waves falls off more slowly with time than the 
corresponding background energy density of the dust driving the FRW expansion. This work 
was subsequently extended by Ressell and Turner |30| who examined the effect of a "dustlike" 
phase during which the scalar field oscillated and decayed on the gravitational radiation 
power-spectrum. The work was then further generalized by Sahni [31|, who repeated these 
calculations for power-law inflation. 

Interest in this subject was re-awakened by the publication of the COBE data ||. A 
number of papers have examined whether the different /—dependence of the scalar and tensor 
contributions to (af) permit one to determine their separate amplitudes. Typically these 
compare the scalar contributions expected from a period of quasi-exponential (slow-roll) 
expansion (which inflates any early perturbations to well beyond today's Hubble radius) 
to the tensor perturbations. These include work by Souradeep and Sahni |p2| , Liddle and 
33H, Davis et. al. |53J, Salopek |55[] , Lucchin, Matarrese and Mollerach pEfl , Dolgov 

and Crittenden et. al. jff5l39|. Another possibility is that one 



Lyth pm , Davis et. al. 
and Silk |37|| , Turner 
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may distinguish the scalar and tensor contributions to the multipole moments by examining 
the polarization of the CBR. This has been examined by Harari and Zaldarriaga j|(J, by 
Crittenden, Davis and Steinhardt [09(1 , and by Ng and Ng pT |. 

Krauss and White have used statistical methods and the COBE data to put tighter 
constraints on the energy- density during an inflationary epoch. Further details of a Monte- 
Carlo simulation were given by White |13[] who also presented a concise derivation of the 
formula for the (af) due to tensor perturbations, and a table of the first ten (of). The effects 
of cosmic variance on the ability to distinguish the scalar and tensor perturbations and the 
slope of the power-spectrum was also considered by White, Krauss and Silk ||4] . 

A related analysis has been performed by Bond et. al. and by Crittenden et. al. [T3j 



who investigate how well one can measure a number of important cosmological parameters 
from the collection of anisotropy observations. It turns out that the first few multipole 
moments are sensitive to very-long wavelength modes which probe well outside our current 
Hubble radius. Stevens, Scott and Silk |46| and Starobinsky |26j have used these measure- 
ments to put new lower limits on the "circumference" of the universe, in the case where it 
has toroidal spatial topology. Such analysis may also be possible in the spatially open case, 
where the Sachs- Wolfe effect has been studied by Ratra and Peebles . 

Grishchuk has also examined the multipole moments arising from gravitational wave 
perturbations |I?|fEJ, adapting the terminology and techniques of quantum optics to carry 
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out the analysis. Grischchuk stresses the importance of the phase correlations between the 
modes of the metric perturbations; we agree with this conclusion but do not use Grischchuk's 
"squeezed state" representation of the field operator. It is possible to obtain identical results 



using only the standard formalism of curved-space quantum field theory developed in [21 



Our conclusion is that the standard ansatz for the so-called "Gaussian spectrum of initial 
perturbations" only gives reliable results for small values of Z; for the higher moments the 
phase relationship between the positive- and negative-frequency components of the wave- 
functions do affect the multipole moments. 

Deviations from Gaussian behavior may in principle be observed through the three-point 
angular correlation function. This was first calculated by Falk, Rangarajan and Srednicki 
SOfl ; the implications of these results and further analysis have been carried out by Luo and 



Schramm [|5T| and Srednicki [[52 

The physical processes giving rise to the CBR temperature fluctuations may be under- 
stood (and explained) in several ways. We repeat the interpretation given by Allen |2T)f . 
which also sheds light on our technical methods. The period of exponential expansion is an 
unstable one, from the global point of view. During this expansion, perturbations of the 
spatial geometry tend to freeze in dimensionless amplitude, so that when viewed globally the 
spatial sections become more and more distorted. However, another consequence of the rapid 
expansion is that locally, any observer can only see (within her Hubble radius) a smaller 
and smaller region of this spatial section. Hence from the observer's local point of view, the 
spacetime is getting closer and closer to a perturbation-free de Sitter spacetime. One con- 
sequence of this global instability/local stability is that gravitational perturbations which 
are of local origin (for example, due to thermal fluctuations) are very rapidly redshifted in 
wavelength and amplitude. At late times, after sufficient inflation, these perturbations are 
no longer visible to an observer; the only perturbations which remain visible are those of 
quantum origin (the zero-point fluctuations associated with the uncertainty principle) be- 
cause these fluctuations extend up to arbitrarily high frequency and can not be redshifted 
away. (In similar fashion, the quanta radiated by an evaporating black hole at late times 
are due to quantum zero-point fluctuations at very high frequency close to the horizon.) 
Hence, to determine the gravitational perturbations present at late times, we assume that 
the initial state of the universe was the vacuum state appropriate to de Sitter space, con- 
taining only the quantum fluctuations and no additional excitations. For this reason, one 
can do a calculation based entirely on "first principles"; the amplitude of the primordial 
fluctuations follows directly from the canonical commutation relations obeyed by the lin- 
earized gravitational field, or in physical terms, directly from the uncertainty principle. The 
"particle production" in this case is the production of pairs of gravitons, whose collective 
effects (since the occupation numbers are large, and they are bosons) appear as classical 
gravitational radiation. Thus, we determine the expected value of the correlation function 
( |1 . ip by finding the expectation value of ^(u a )^-{v h ) in the de Sitter vacuum state. 

The published calculations have two shortcomings which are addressed in the present 
work. The first is pedagogic. The calculations which have been published are all rather 
sketchy; to reproduce the results requires many pages of calculation which are not given in 
full but are left as an exercise to the reader. We believe that our method of performing 
this calculation is new; it is short enough and elegant enough so that all of the details 
can be shown explicitly. The second advantage is quantitative. The previously published 
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calculations are made using a "long wavelength" approximation to the mode functions, which 
is accurate for determining the effects of the lowest multipoles, but inaccurate for the higher 
ones. It is not obvious from the published work how to remove this approximation to obtain 
more accurate results; in the present work we give exact expressions for the correlation 
function. The shortcomings of the standard "long wavelength" approximations have been 
pointed out in the recent work of Turner, White and Lidsey f53f , who use numerical methods 
to integrate the wave equation for the mode functions and who obtain results that appear 
identical to our exact formula. 

This paper is organized as follows. After a few notes on notation, section [TT] begins with 
the classic formula for the Sachs- Wolfe effect in spatially flat FRW cosmological models. 
This is used to derive an expression for the correlation function C^j) due to the gravita- 
tional radiation, under the assumption that the initial state of the universe is a vacuum 
state with only zero-point perturbations. In section |T| we derive from first principles the 
normalization condition on the graviton wave functions. Section [TV] describes a simple infla- 



tionary cosmological model, also used in p9[ . In this model, the universe "begins" with an 
infinite period of inflation, then makes an instantaneous transition to a radiation-dominated 
stage, and then later makes another instantaneous transition to a matter-dominated stage. 
We then find the normalized graviton wave functions appropriate to that model (and the 
corresponding Bogolubov coefficients). We also show how the standard results appear as a 
low-frequency approximation to the exact expressions. Section [V| is an attempt to obtain 
a closed form for C(7); this attempt does not succeed but some progress is made. Section 
V| compares the results of our exact expression for the expansion coefficients (of) with the 
more standard results, and includes a discussion of some recent literature on the subject. 
Because the high-frequency modes affect the temperature perturbations on small angular 
scales, the exact (af) agree with those given by the standard approximations for small /, and 
are different for large I. Finally a pair of appendices show an alternative derivation of the 
formulae contained in section |TJ, and contain a brief description of the numerical techniques 



used in section VI 



Throughout this paper, we use units where the speed of light c = 1. However for clarity 
we have retained Newton's gravitational constant G and Planck's constant h explicitly. 



II. THE SACHS- WOLFE EFFECT 
AND THE ANGULAR CORRELATION FUNCTION 

A. Notes on Notation 

We begin with a few notes on notation. The vectors and tensors in this section are purely 
spatial; they have no time components, although they may be time-dependent functions. In 
a spatially-flat FRW model, the spatial geometry is flat Euclidean space. Since the tensors 
and vectors are spatial we raise and lower tangent space indices with the spatial part of the 
conformal metric, which is just the Euclidean metric of 3? 3 . In Cartesian coordinates, this is 

<U = diag(l,l,l). (2.1) 

We denote spatial vectors by k a , v a , or u a , and spatial tensors by h ab , e ah , or P ab . The Latin 
indices a,b, . . . , f run from 1 to 3. Associated with any spatial vector is its magnitude, 
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denoted by the vector symbol without a tangent space index. For example, the magnitude 
of the vector k c is denoted k, where 



k = yfk"k a = sJ5 ab k a k\ (2.2) 

A special notation is used for spatial vectors with unit magnitude. The unit spatial vector 
k a is defined by 

k a 

k a = (2.3) 

so that k a k a = 1. Thus one may decompose any spatial vector u a into a magnitude and a 
unit vector, and express it as 

u a = uu a . (2.4) 

We use this notation throughout this section. 

Often we need to integrate over all possible magnitudes and orientations of a spatial 
vector. To integrate over k a , we write 

p roc /*7r p2tv roo r 

Jd 3 k = J dkk 2 J d8 k sin8 k J d<p k = J dk k 2 J dQ k , (2.5) 

denoting the polar angle associated with k c by 8 k , and the azimuthal angle by <p k . In a 
similar way we will denote a function of the polar and azimuthal angles (such as a spherical 
harmomic function) as 

Y lm (0k,<Pk)=Y lm (k c ). (2.6) 
For example, the orthonormality condition for the spherical harmonics is 

J dn % Y lm {k c )Y; q {k c ) = 5i p 5 mq . (2.7) 

Note that we never integrate over the polar and azimuthal angles separately. 
Hilbert space operators are denoted by an overbar, for example 

<#> = m, (2.8) 

and a f denotes the adjoint operator. A * denotes complex conjugation. 



B. The Sachs- Wolfe Effect 

In a perfectly isotropic universe the CBR would have the same temperature in all di- 
rections on the celestial sphere. If, however, the cosmological metric is perturbed away 
from isotropy, the temperature observed today fluctuates over the celestial sphere, even if 
the last-scattering surface had uniform temperature. The Sachs- Wolfe formula expresses 
the temperature fluctuations in terms of the metric perturbations, to lowest order in the 
perturbation. 



S 



Consider a spatially-flat FRW universe perturbed away from isotropy. The perturbed 
metric in comoving gauge is 



ds 2 = a 2 (7])[-drj 2 + (S ab + h ab (r), x c ))dx a dx% (2.9) 

where n is the conformal time, and a(rj) is the scale factor. Imagine a single photon emitted 
at conformal time rj = rj e and observed at conformal time rj = r/ b s . As the photon propagates 
through the spacetime the metric perturbations a 2 (r])h a b perturb the null path of the photon. 
One may choose spatial coordinates so that (to zero order in h ab ) the photon follows a radial 
path in space on its way to the observer, who is located at the origin x a = 0. Furthermore, 
one may parametrize the path of the photon by A, so that the spatial path of the photon is 

x a (X) = D(X)u a , (2.10) 

where 

D{X) = (rfabs —r] e — X), (2.11) 
u a is a unit vector pointing radially out from the origin, and A varies from A e to A Q b s with 

A c = 0, (2.12) 

Aobs = Vobs - Ve- (2-13) 



Sachs and Wolfe have shown that to first order in h ab the observed redshift of the photon is 
given by 



1 + Z 



a(r], 



obs I 



a(r] c 



1 



2 Jx 



u a u b 



d_ 

dr] 



h ab (r],D(X)u c 



dX). 



?7=?7e+A 



(2.14) 



This equation is equivalent to (39) in reference for the specialized case of gravitational- 
wave perturbations. 

The CBR is an ensemble of many photons which were last scattered at conformal time 
rj = r] e by the primordial plasma of ionized hydrogen and electrons. Using fl2.14|) one obtains 
the temperature fluctuation ST of the CBR measured at the point on the celestial sphere 
pointed to by the unit vector u c : 



ST 
~T 



1 

2J\ 



A,-, 



u a u b 



d_ 

dr] 



h ab (rj,D(X)u c 



dX. 



(2.15) 



This formula embodies the Sachs- Wolfe effect, and is equivalent to (42) in reference for 
the special case of gravitational-wave perturbations. 



C. The Metric Perturbation h ab 

As noted in the introduction, we examine the transverse, traceless, tensor part of the 
metric perturbation in models of the universe that pass through an early inflationary stage. 
The period of exponential inflation is unstable, and as a result of the rapid expansion, 
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perturbations of the spatial geometry freeze in dimensionless amplitude. From any observer's 
local point of view the spacetime quickly approaches a perturbation-free de Sitter spacetime. 
At late times perturbations of local origin are extremely redshifted in both wavelength and 
amplitude, leaving only perturbations of quantum origin (zero-point fluctuations) as the 
significant contribution to the tensor part of the metric perturbations. For this reason we 
assume that the initial state of the universe is the de Sitter space vacuum state containing 
only quantum fluctuations. 

Since the significant tensor perturbations are quantum in origin, we replace the classical 
metric perturbation h ab in ( p.!5| ) by the Hilbert space operator h ab appropriate for the 
linearized theory of gravity. The plane wave expansion of h ab is 

h ab (v,x c ) = J d?k (e ikdx *[ eab (k c )Mv,k c )a R (k c ) +e* ab (k^)Mv,k e )a L (^)] 

+e-^[e: b (A; c )0^(r / ,A ; c )4(A; c ) + e ab (k^U V , k c )a{{k c )]). (2.16) 

Here dji(k c ) and a^(/c c ) (their Hermitian conjugates) are annihilation (creation) operators 
that destroy (create) a right or left circularly polarized graviton. These operators obey the 
commutation relations 

[a L (k a ), a[(k' a )] = [a R (k a ), a R (k' a )] = 5 3 (k a - k' a ), (2.17) 

with all other commutators vanishing. The graviton mode functions for the left and right 
polarizations are 4>L(rj,k c ) and 4>n(ri,k c ) respectively. If the spacetime is isotropic and ho- 
mogeneous, and therefore does not single out any preferred directions, one may choose a 
particle basis so that the mode functions depend on the magnitude k only. One may also 
choose a particle basis that does not distinguish between the two possible spatial orienta- 
tions, so that the left- and right-handed gravitons have the same mode functions. One then 
has 



Mv,k c )=Mv,k c )=^V,k). (2.18) 



This mode function (f)(r] } k) obeys the massless Klein-Gordon equation [pT 



a(rj) 

where a(rj) is the cosmic scale factor, and 



+ 2 ^U + fc 2 = O, (2.19) 



d_ 

dr] 



(2.20) 



If one demands that h ab obey canonical commutation relations, the commutation relations 
( p. 17] ) imply that the mode function satisfy normalization conditions. The normalization 
condition is defined in ( |3.19|) . 

The tensors e a b(k c ) and e* ab {k c ) in the expansion Q2.16Q are the polarization tensors for a 
circularly polarized basis. We first define the so called plus (+) and cross (x) polarizations. 
Consider a mode or wave propagating in the k c direction. One may define two unit length 
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vectors rh c and n c orthoganol to k , and orthoganol to each other, so that the set (k c , rh c , h c ) 
is a right-handed triad with 

k a m a = k a h a = m a h a = 0. (2.21) 

In terms of these unit vectors the plus and cross polarizations are defined as 

e ( aV(k c ) = m a (k c )m b (k c ) ~ n a (k c )h b (k c ) (2.22) 
e%\k e ) = m a (k c )h b (k c ) + n a (k c )m b (k c ). (2.23) 

The plus and cross polarization tensors together form a complete basis for the tensor (spin 2) 
perturbations |J. Note that both the plus and cross polarizations are transverse, traceless, 
and symmetric: 

e£\k c )k a = e^\(k c ) = e[ij(A; c ) = 0. (2.24) 

One may define the circular polarization tensor e ab (k c ) in terms of the plus and cross polar- 
izations as 

e ab (k c ) = l=[ei + b \k c )+zei x b \k% (2.25) 
= -^[m a (k c ) +tn a (k c )}[m b (k c ) +in b (k c )\. (2.26) 

The polarization tensor e* ab (k c ) is just the complex conjugate of the polarization tensor 
( |2.25| ). The tensors e ab (k c ) and e* ab {k c ) also form a complete basis for the tensor (spin 2) 
perturbation. 

The vectors m c (k c ) and n c (k c ) are not unique. Any two unit vectors that satisfy ( |2.21D 
may be used to define the polarization tensors. Any other right-handed triad of vectors such 
as (k c , rh' c , n' c ), however, can be obtained by rotating the triad (k c ,rh c ,h c ) through an angle 
<p about k c . Under this rotation, 

e' ab (k c )=e- 2 ^e ab (k c ). (2.27) 

This shows that gravitons are a spin 2 field, since the spin (or more precisely, the helicity) 
of a field is defined as the number of times the phase of the field changes by 27r, when the 
coordinate system is rotated once around the momentum vector of the field. 

The polarization tensors are closely related to the tensor that projects onto a sphere of 
radius k, at a point k c . We define the projection tensor P ab {k c ) by 

P ab (k c ) = 5 ab - k a k, (2.28) 

so that 

P ab k a = P ab k b = 0, and P ab P b c = P ac . (2.29) 

This tensor projects onto the two-surface orthoganol to k c , which is just the two-sphere 
of radius k. To relate the projection tensor P ab to the polarization tensors, consider the 
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Euclidean metric 6 ab on 3? 3 . One may express 6 ab using the three unit vectors k c ,rh c , and 



Sab = k a h + rh a rh b + n a n b . (2.30) 
Using ( |2.28| ) and ( |2.30| ) one may write the projection tensor as 

P ab (k c ) = m a (k c )m b (k c ) + n a {k c )h b {k c ) . (2.31) 
From ( |2.26| ) and ( |2.31D one may quickly verify the identity 

e ab (k e )e: d (k e ) + e: b (k e )e cd (k e ) = P ac (k e ) P bd (k e ) + P ad (k e )P bc (k e ) - P ab (k e )P cd (k e ). (2.32) 
Later we use this identity to find an elegant expression for the angular correlation function. 



D. The Two-Sphere of Radius k 

Besides being the projection tensor onto the two-sphere of radius k, P ab is the natural 
metric induced on this two-surface by the flat metric on 3? 3 . Since the two-sphere is a 
maximally symmetric two-manifold, one may immediately write the Riemann tensor on this 
two- surface as 



Pabcd ,2^ >a [ c ^ > d]b' 



[2.33) 



The factor of k 2 appears because the two-sphere has radius k. We denote the covariant 
derivative on this surface by V , and define the Laplacian □ on this surface by 



P ah V a V b 



v b v. 



□. 



(2.34) 



The spherical harmonics are eigenfunctions of this Laplacian, and obey the eigenfunction 
equation 



k 2 



(2.35) 



Again the factor of k 2 appears because the the two-sphere has radius k. 

Using the definition of the Riemann tensor, the identity ( |2.33| ), and the eigenfunction 



equation ( |2.35| ), we can derive a useful identity: 

7 y* 

b 1 Ira 



(Vv a v 6 - v a v b v b + v a v b v b )y;^ 

R ba bc V%* m + V a aY; m 
—\7 a Y,* +V a DY? 

^,2 ™ ' m 



-1(1 + 1) + 1 



k 2 



(2.36) 



This formula will prove useful in our derivation of the angular correlation function. 
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E. The Angular Correlation Function 



1. The Sachs-Wolfe Operator 



The Sachs- Wolfe formula ( 2.15 ) is a result from classical general relativity, giving the 
temperature fluctuations of the CBR over the celestial sphere as a function of metric pertur- 
bations. As noted above, however, in inflationary models the surviving metric perturbations 
are quantum in origin; without further justification we replace the classical metric perturba- 
tion h ab in the standard Sachs- Wolfe formula ( [2.15| ) by the quantum field operator h ab . The 
temperature fluctuation at a point on the celestial sphere is now a Hilbert space operator, 
given by 



5T 
~T 



1 Mobs 

2 A c 



u a u b -^h ab (v,D(\)u c ) 



dX. 



)?=)?e+A 



(2.37) 



We will refer to (|2.37|) as the Sachs- Wolfe operator. 

Since the Sachs- Wolfe operator is parametrized by coordinates on the celestial two-sphere, 
it is natural to decompose it into an expansion of (normalized) spherical harmonics on the 
two-sphere. Using the orthogonality of the spherical harmonics, one can write the Sachs- 
Wolfe operator as 



5 -^(u a ) = T,Ci m Y lm (u a ) 

Im 



(2.38) 



where the expansion coefficient operator Ci m is 



a 



i r x - 



Im 



2 7a 



dX I dQ^u^Y^iu ) 



d_ 

dr] 



h ab (r,,D(X)u c ) 



7]=7] e +X 



[2.39) 



For the metric perturbation operator h ab we use ( p,16| ). Since the first derivative of the 
metric perturbation, not the perturbation itself, appears in the expression for the expansion 
coefficient Ci m , it is useful to define the dimensionless function 



F(X, k) = k 1/2 



d_ 

dr] 



(f){r],k) 



(2.40) 



Then from ( |2.16| ) and ( |2.39| ) we obtain for the expansion coefficient operator 

_ A ba 



a 



Im 



dX I dVtt 



d 3 k 

Ac 

x|e^W fc ^F(A,A;)[e a6 (A: c )afl(^) + ei(A: c )o L (A: c )] 
+e- w ^ kd ^F*(X, k)[e: b (kTaUk c ) + e ab (k c )a[(k c )] 



(2.41) 



We use this expansion of the Sachs- Wolfe operator to examine the angular correlation func- 
tion. 
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2. Angular Correlation Function C(v a ,u a ) 



The quantity of interest is the angular correlation function (|1.1|) . Since the temperature 
fluctuations are now represented by a Hilbert space operator, the angular correlation function 
is a matrix element: 



C{v a ,u a ) 



. 



^2.42) 



Here the quantum state |0) is the initial quantum state of the universe, which we have taken 
to be the de Sitter space vacuum state for reasons discussed both in the introduction and 
in section 111 C. 



Based on the isotropy of the FRW model and of the state |0), one expects the correlation 
function to be rotationally invariant; i.e. to depend only on the angle 7 where cos 7 = v c u c . 
Using the expansion ( p. 38 ), one may express the angular correlation function in the form 



im pq 



(2.43) 



Since one expects the correlation function to be rotationally invariant, one ought to be able 
to write the matrix element (0\C^ q Ci m \0) as 



<0|Cj.Ci m |0) = {a?)6 lp 6.. 



mqi 



(2.44) 



and then use ( |2.41j ) for Ci m and solve for (af). In appendix |A| we make this assumption, 
and obtain (of) somewhat more directly. 

For now, however, we show by direct calculation that the correlation function is rota- 
tionally invariant. Using fl2.41|) the matrix element (0|(3 q Ci m \0) is 



(0\& pq C lm \0) = - 



dX' / dX 



d 3 k' r d 3 k 



k' 1 / 2 J k 1 / 



—F(X',k>)F*(X,k) 



e afc (A:' e )e:,(F)(0|a fi (^)4(A; e )|0) + e: fe (A;' e )e cd (F)(0|a L (A;' e )4(A; e )|0) 

x J dtlt J ffii^^^^^l^^r^e-^^W^- ^''^'. (2.45) 

One may immediately evaluate the two matrix elements on the right-hand side using the 
commutation relations ( |2.17| ) for the creation and annihilation operators. Both matrix ele- 
ments yield the Dirac delta function 5 3 (k c — k' c ). Using the identity ( |2.32j ) for the polarization 
tensors, one finds 

(0\Cl q C lm \0) = \j*° hS dX' £° bS dxj™ dkkF(X',k)F*(X,k)A lmpq (k,D(X),D(X')), (2.46) 
where 

A lmpq (k,r,r') = J da k {[P ac (kk e )PUkk e ) + P ad {kk e )P bc {kk e ) - P ab (kk £ )P cd (kk e )} 

xr4 q} (r'kk e )rtL } (rkk e )\, (2-47) 
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and 

*P{L}(k c ) = J dfi<i y im (u c )«V fc ^. (2.48) 

The braces in the equation above are to remind the reader that I and m are not tangent space 
indices. We show in the next section that Ai mpq (k,r,r') is proportional to the Kronecker 
deltas Si p S mq and is independent of m, so that the correlation function is indeed rotationally 
invariant. 



3. A Closed Form Expression for A[ mp g(k, r, r') 

The function ^i m y( k c ) can be expressed in a way which allows one to exploit the pro- 
jection tensors in ( 2.47|) . Note from ( |2.48| ) that 



^{L } (k c ) = -V a V b J cM* Y lm (u c )e tkd ^ } (2.49) 
where the derivative V a in Cartesian coordinates is 

V.^. (2-50) 

The plane wave e lkdud can be expanded as an infinite sum of spherical Bessel functions ji(k) 
and spherical harmonics plf so that 



P r P 

^ m} (k c ) = -V"V b / <Kl & Y lm {u c ) ^J2i p j P (k) E Y^Y^V 

p=0 q=— p 

Using the orthonormality of the spherical harmonics one obtains 

^{l}(k c ) = -Am l V a V b 3l (k)Y lm (k c )- (2.52) 

Note the dependence of the right-hand side on the vector k c ; the spherical Bessel function 
depends only on the magnitude k, and the spherical harmonic depends only on the polar 
and azimuthal angles. With this form for ip?i m }{k c ) ( |2.47| ) becomes 



(2.51) 



Almpq{k } 7", T ) 



l6n 2 iP( 



-i 



x J dn k l [ [p ac (kk e )p bd (kk e ) + p ad (kk e )p bc (kk e ) - p ab (kk e )p cd {kk e )} 

x[V c V dJp (A;r')r M (A ; e )][V a V^(^)^(^)]}- (2-53) 

We can now use the projection operators to make the final integration almost trivial. 

Consider how the projection tensor P a b acts on the gradient V a f{k c ). The gradient 
in general has components both parallel and orthoganol to k c . When contracted with the 
gradient the projection tensor annihilates the components parallel to k c . The remaining 
components of the gradient lie entirely on the two-sphere of radius k, so 
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P\V a f(k c ) = V b f(k c 



(2.54) 



where V a is the same derivative on the two-sphere defined in section |11 D , Using ( |2.54| ), and 
noting that the spherical Bessel functions depend only on the magnitude k, and are constant 
on the two-sphere of radius k, one obtains 



16ir 2 iP(-i) 1 



iy 2,y/2 



j P (kr')ji(kr) 



x / dn~ k [2(v a v b Y pq )(v a v b Y;j - (□y M )(n^)], 



(2.55) 



where □ is the same Laplacian on the two-sphere of radius k defined in (|2.34|) . The integrand 
is just derivatives on the two-sphere of spherical harmonics, which, as discussed in section 
H D| , are eigenf unctions of the Laplacian □. 

The first integral on the right-hand side above can be integrated by parts, and the second 
by inspection. To help us evaluate the integrals, we write 



Aimpq(k, r, r ) 



16tT 2 ^( 



ly 2 ^"/2 



where 



and 



= / ^(□F w )(DF i 



Jp (kr%(kr)[2Q\ 2 2 pq (k) - Q\2 pq (k^ 
l 2 (l + l) 2 



lm) 



k A 



Slr,5 



lp"mq i 



Qimpq(k) = I ^(v a v fe y M )(v a v^ 



lm) 



Integrating (|2.58|) by parts once, we find 



= - / dn % (v b Y vq )(uv% 



b\r* 



lm)i 



since the two-sphere has no boundary. With the identity (|2.36|) we have 



-/(/ + !) + ! 



A- 2 



dn- k (v b Y pq )(v b Y;j. 



(2.56) 



(2.57) 



(2.58) 



(2.59) 



(2.60) 



One may again integrate by parts and use the eigenfunction equation for the spherical 
harmonics to obtain 



Q\mpq{^) 



-1(1 + I) + 1 

k 2 



-1(1 + 1) 



A- 2 



$lp&mqi 



-[l 2 (l + l) 2 -1(1 + 1)} 8 lp 8, 



mq- 



Substituting ( ^57|) and (|2T62|) into ( ^56|) one has 

16tt 2 



Ai mp q(k, r, r ) 



ji(kr%(kr) (I - 1)1(1 + + 2) 5 lp 8 mq . 



(2.61) 
(2.62) 

(2.63) 



As previously indicated, the matrix element (^C^Cim^S) and the correlation function 
C(v c , u c ) are indeed rotationally invariant. Also note that ( |2.63| ) vanishes for I = and 
I = 1. 
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4- The Angular Correlation Function 



Using the above form of Ai mpq (k, r, r') one may derive a simple expression for the angular 
correlation function C(v c , u c ), and show directly that it depends only on the angle 7 between 
v c and u c . With (|2.63|) and ( [2.46| ) one has 



(0\Cl q C lm \0) = ^^-^5 lp 5 mq 



J ~k J J ' (A ' fc) fc2^(A)^(A0 ' (2 - 64) 

where we have written the fourth-order polynomial in I appearing in ( 2.63| ) as the ratio of 
two factorials. Noting the symmetry of the right-hand side, and recalling the definitions of 
D(\), A e , and A Q bs, we define 



rv 0bs -r, c ^ ^ ^ Jl{k{7] obs -T] e - A )) 

k(r]obs ~ Ve - A) 

and write the matrix element as 



!,(*) = /_ d\ F(X, k) J T,: ' 0bS „ e \7 , (2-65) 



P9 

where 



(0|Cj-CyO) = (a z 2 )M m<?) (2-66) 



Substituting this expression into ( p.43| ) we obtain 

C(v c , u c ) = Ei^YirnmY^n (2.68) 

I in 

where we have used the Kronecker deltas to eliminate two of the sums. Making use of the 
addition theorem for spherical harmonics (see equation 3.62 in reference |j54]| ), the correlation 
function is 

00 (21 + 1) 

C(v c ,u c ) = C(j) = £^-^(a 2 >PKcos 7 ), (2.69) 
1=0 47r 

where 

cos7 = m c £ c . (2.70) 

As promised, the angular correlation function depends only on the angle 7 between any two 
points on the celestial sphere. Also note that the I = and the I = 1 terms in the expansion 
vanish exactly. 

This form of the correlation function is very general. The only dependence of the cor- 
relation function on the details of any cosmological model is through the graviton mode 
function (or more precisely, its first derivative), which appears as F(X,k) in the definition 
of h{k) (|2.65 ). A similar result, which is as general as ( |2.69| ), is given by Grishchuk (see 
equation 4 in the second paper of reference |E|). 
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III. GRAVITON MODE FUNCTION NORMALIZATION 



If one demands that the metric perturbation field operator h ab obey canonical commu- 
tation relations, the commutation relations ( |2.17| ) for the graviton creation and annihilation 



operators imply that the graviton mode function satisfy a normalization condition. Impos- 
ing canonical commutation relations on the tensor field h a b, however, is subtle because as 



noted by Ford and Parker [21], the canonical commutation relation that h a b obeys may be 



inconsistent with the gauge conditions on h a b. For this reason, we follow reference and 
impose canonical commutation relations on the two independent scalar degrees of freedom 
in hoi,. 

The two independent degrees of freedom in the metric perturbation field can be isolated 
by constructing two scalar field operators from h a b. Recall from ( 2.16Q that the plane wave 
expansion of h a b is 



hab(v,x c ) = J ^kKbi^.x^k ), (3.1) 

where 

h ab ( V ,x c ,k c ) = (e ikdx *{e ab (k c )Mv,k c )a R (k c ) + e* ab {k c )<f> L (v, k c )a L (k c )} 

+e-^[e: b (A; c )0^r / ,^)4(^)+ e ab (k c )<f>* L ( V , k^aW)]). (3.2) 
We define the scalar field operator 

h+(r), x c ) = J d 3 k h ab {r], x c , k c )e i+)ab {k c ). (3.3) 
Contracting the integrand using ( |2.21| - |2~2~6| ) one obtains 



h + ( V ,x c ) = V2j al 3 k |e* fc Xd (p(i], k)\d R {k c ) + a^(A; c )] 

+e- ikdx ^*( V ,k)[a R (k c ) +a L (fc c )] t }. (3.4) 

A second scalar field operator h x (r],x c ) is defined by replacing the plus signs (+) in fl3~3[) by 
crosses (x), which has the effect of replacing a^+a^ by ia R —ia L in (|3.4|). Together the scalar 
field operators h+(i], x c ) and h x (r/, x c ) possess the same two degrees of freedom as the metric 
perturbation operator h a b pl| . Since h + (r],x c ) and h x (i],x c ) are both scalar field operators, 
they obey well known canonical commutation relations for scalar fields. Because our particle 
basis does not distinguish between the two polarizations, we only need to consider one of the 
two scalar fields, since both lead to the same normalization condition for the mode function. 
The scalar field operator h + (i],x c ) obeys the canonical commutation relation 

[h+(r], x c ), 7r+(77, x' c )} = ih5 3 {x c - x' c ), (3.5) 

where the field operator 7t + (r],x c ) is canonically conjugate to h + (rj,x c ), and is defined by 

vr + ^ I*. (3.6) 
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Here C is the Lagrangian density of the perturbed spatially-flat FRW spacetime. To impose 



the commutation relation (3.5) and find the normalization condition for the graviton mode 
function we need to find the Lagrangian density in terms of h + and h x . 

The Lagrangian density is obtained by expanding the gravitational plus matter action 
to second order in the metric perturbation h a b- The action for a FRW spacetime is 

S = |rf 4 x v ^{ I ^_ + i[(p + PK^ + (p + 3P)]}, (3.7) 

where is the metric for the FRW spacetime, R is the Ricci scalar, p is the energy density, 
P is the pressure , and is the four velocity of the cosmological fluid. Varying the action 
with respect to g^ u leads to the Einstein equation for a FRW model: 

= R P u + ^Rg^u = SttGT^ = 8vrG[(p + P)u^u u + Pg^\. (3.8) 
If the FRW spacetime is perturbed so that 

fa = ° V + ( 3 - 9 ) 

where °g IMU is the unperturbed or background FRW metric, but the pressure P and the 
energy density p are not perturbed, then one finds that the second-order variation in the 
action is |2T] 



5 2 S = J ^V^^{(°V^)(°V M 7^) + SnG(P - p)^^ 

+2°/iV7 M V + 2°iW7 M<7 7^}- (3-10) 

The superscript, for example in °V M , refers to the background spacetime. One should note 
that ( p.lO| ) is obtained by making a specific choice of gauge (transverse, traceless) pip . Also 
note that this is to second order in the perturbation j^, since the first order part 5S vanishes 



because the background FRW spacetime satisfies ( |3.8| ) 



Equation ( [3.10| ) is very general and true for any "small" perturbation 7^ (that satisfies 
the gauge conditions) away from a FRW spacetime with metric °g liV . For our purposes, the 
perturbation 7^ is simply a 2 {rj)h a i,, and the background FRW spacetime is spatially-flat. 
With a little calculation one can show that for the spatially-flat FRW spacetime perturbed 
by tensor perturbations 

S 2 S = J - h ab h ab + {d a h bc Wh hc )Y (3.11) 

To calculate the momentum 7f + conjugate to h + , one must express the action in terms of h + 
and h x . With a little calculation, and using (|2.16| ), one can write the action as 



S 2 S = J d 4 x^^{-(h + + h x ) + (d a h + )(d a h + ) + (d a h x )(d a h x )}, (3.12) 
so that the Lagrangian density is 
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C = Wg\ { ~^ + + ^ + (W^) + (^x)(9^ x )}. (3.13) 

This is just the Lagrangian for a pair of massless scalar fields minimally coupled to the 
background spacetime. Using (3J3) and ( |3.13|) the momentum is 



= -iis^. (3 - 14) 

Then from ( |3.5| ) the canonical commutation relation for the scalar field operator h + is 

_. r) 3 (r a T ,a ) 

[h+(r], x a ), h+(r], x' a )] = -Mnih G 1 -. . ; . (3.15) 

0/(77) 

Note that this is an equal-time commutation relation. 

Using the commutation relation above and the explicit form for the scalar field h + , one 
can derive the normalization condition for the graviton mode function. Using (ft.4|) and 
( p,17| ) one finds 



[h+fa x a ), h + (rj, x' a )] = Aj d 3 k{(j)(r}, k)4>*(r}, k)^^ - 0*(r/, k)^, k)e- ik ^ x *- x '^} . 

(3.16) 

Since we assume that the mode function 0(r7, k) depends only on the magnitude k, one can 
write 

[h+{n, x a ), h+(r], x' a )} = 4 / d 3 ke ikaiXa - x '^{(f){r], typfa k) - <j)*{r], k)j>(ri, k)}. (3.17) 



Since the delta function in (|3.15 ) can be expressed as a plane-wave expansion 



d 3 ke ik"( Xa -x> a ) = (2n) 3 5 3 (x c - x' c ), (3.18) 
( |3.15| ) and ( |3.17| ) imply the mode function normalization condition 

{<p( v , k)4>*(i], k) - <f>*{ v , k)fa, k)} = -^- v (3.19) 

This identity determines the normalization of the graviton mode function, up to an (irrel- 
evant) overall phase, and would be equivalent to equation (3.3) of reference |2l] if not for 
a typo [53] . The main consequence is that fundamental physical principles (the uncertainty 
principle) completely determine the amplitude of the contribution to the angular correlation 
function arising from gravitational radiation. 



IV. INFLATIONARY COSMOLOGICAL MODEL 

A. Graviton Mode Function 

The cosmological model we examine "begins" with an infinite inflationary phase, followed 
by radiation- and then matter-dominated phases. We assume that the inflationary phase 
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evolves into de Sitter spacetime. The mechanism by which the universe arrives at the de 
Sitter spacetime is not important, since the period of rapid expansion during the de Sitter 
phase effectively erases the initial conditions. At the end of the de Sitter phase the universe 
undergoes an instantaneous phase transition to a radiation-dominated FRW phase. At the 
end of the radiation phase the universe again undergoes an instantaneous phase transition, 
and evolves as a matter-dominated FRW spacetime until the present. 

If the initial de Sitter phase is sufficiently long, the spatial geometry becomes flat, and 
one may assume that the universe is spatially-flat for all three epochs. The metric for the 
spacetime is then given by (|2.9|), with scale factor 



a{r]) 



(2 — -^-) 1 a{rji) — oo < 7] < r]i deSitter, 

^-a(?7i) rji < r] < r] 2 radiation, (4.1) 

„ K 1 + m^m a ^ ^' 2<r l matter, 

where r]i and r]2 are constants. The redshift at the end of the de Sitter phase Z en d and the 
redshift at the time of radiation-matter equality Z e quai are defined by 

1 + Z cnd = — — - = , (4.2J 

0(771) 4??i?72 

1 i_ y - a (Vobs) _ (Vobs + mf on 
J- + ^equal = / r" — Z o ' 

a{r] 2 ) Ar] 2 2 

where r/ bs is conformal time today. Typical values for the redshifts (for models "with enough 
inflation" to solve the horizon and flatness problems) are Z en( j ~ 10 27 and Z eqaa \ w 10 4 . We 
assume that last-scattering at conformal time rj e took place after the time of radiation-matter 
equality so that r\ e > r/ 2 - The reshift of the surface of last-scattering Z\ s is 

l + Z ls = ^4=( r] ° hs + r]2 ) 2 . (4.4) 



a(r] c ) V rj c + 7] 2 

A typical value for the redshift of the surface of last-scattering is Z\ s w 1300, although it is 
possible that the hydrogen was re- ionized as recently as redshift Z ls w 100 [|]. 

Note that the scale factor ( |4.1|) and its first derivative are continuous. Because the 
second derivative of the scale factor is not continuous, the scalar curvature of the spacetime 
changes discontinuously at the phase transitions. This instantaneous phase transition is 
a good approximation, except at high frequencies, where it predicts too much graviton 
production |H5| . 



With the scale factor above, one can solve the massless Klein-Gordon equation (p. 19 ) 



for the graviton mode function during each of the three epochs. By making a change of 
dependent, and then independent variable, the Klein-Gordon equation can be cast in the 
form of Bessel's equation, for each of the three phases. The necessary changes of variable, 
and the positive-frequency solutions, are shown in table |. Using the normalization condition 
( P.19| ), and making a convenient choice of phase, one obtains the following positive-frequency 
solutions for the three epochs: 
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^+\ v , k) = -iJ—— k 1 ' 2 (77 - 2 Vl ) 2 hf\k(r] - 2r ll ))e- ik ^\ for the de Sitter phase, (4.5) 

V 3?r Pp 



<j£)(r}, k) = -iJ—— k 1/2 T/i 2 (kf])e ik71 \ for the radiation phase, (4.6) 
&(V, k) = -UJ^-— k 1/2 Vl 2 V2 Hl t^ + ri) for the matter hag (4J) 

\ •)-.' Pp '/ • '1-2 



where 



3 a 2 (Vi) 3 1 



is the (constant) energy density during the de Sitter phase, and 



Pp = ^ (4.9) 



is the Planck energy density. The spherical Hankel functions 54 are defined by 



h\*\z)=j l (z)±iy l (z), (4.10) 

where ji(z) and yi(z) are spherical Bessel functions of the first and second kind. The negative- 
frequency mode functions are the complex conjugates of the positive- frequency mode func- 
tions. The positive- and negative-frequency solutions for each epoch form a complete set of 
solutions to the massless Klein-Gordon equation ( |2. 19|) . 



The choice of a mode function during the initial de Sitter phase 77 < 771 completely 
determines the mode function at all later times. This is because a solution to the Klein- 
Gordon equation (|2.19| ) depends only upon the values of and on a spacelike hypersurface 



(i.e. a surface of fixed rj). To express the solution at later times, after the de Sitter 
phase has ended, it is useful to adopt the Bogolubov coefficient notation. In this notation, 
the solution a ^ later times is expressed as a linear combination a0 < - + - ) + (3$' of the 
natural choices of positive and negative frequency solutions during the subsequent phases 
of expansion. 

If one evolves the positive-frequency mode function during the de Sitter phase 0^ into 
the subsequent radiation phase via (|2.19| ), it is necessary that the mode function and its 
first derivative be continuous across the phase transitions at 77 — 771 and 77 = 7/2. Continuity 
from the de Sitter to the radiation phase is assured if and only if the Bogolubov coefficients 
a ra d and /5 ra d satisfy the conditions 

<f>da {Vi, k) = cwrfaj^i, k) + ^0^(771, k), 

0dt (Vi, k) = a rad 0S(77i, k) + /^^(r/i, k). (4.11) 
Solving this pair of linear equations one finds 



a r ad(^i, k) = -il 1 + 



1 



kr/i 2k 2 r/i 2 

^i,k) = ^- r (4.12) 
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Likewise, if one evolves the positive-frequency mode function during the radiation phase <p^\ 
into the subsequent matter phase, the Bogolubov coefficients a mat and /3 mat must satisfy 



0rad(^2, k) = a mat &t(V2, k) + Anat0mal(r?2, k) . 



(4.13) 



Solving this pair of linear equations we find 

a m at(?72, k) = -i(l + 

Prazt(ri2,k) = — 5— ^ 



2kr/2 8/c 2 r/2 5 

ifc(?jl-3?72) 



3 ifc(?7i+r;2) 



(4.14) 



Since the mode functions are normalized by ( p. 19 ), the Bogolubov coefficients obey the 
(easily verified) relation 



l^radl 2 — |A-ad| 



mat 



1. 



(4.15) 



The Bogolubov coefficients above agree with reference |29[, up to an irrelevant phase. 

As stated earlier, the choice of a mode function during the de Sitter phase completely 
determines the mode function at all later times. We choose the mode function for the de 
Sitter phase to be the positive- frequency de Sitter solution ( |4.5|) : 

(f)(r],k) = (j)fa\r},k), for — 00 < 77 < rji. (4.16) 

This is the unique solution corresponding to a de Sitter-invariant vacuum state with the 
same (Hadamard) short distance behavior as one would find in Minkowski space [29]. Having 
calculated the Bogolubov coefficients, one may now determine the way in which the positive- 
frequency mode function ( |4.16| ) evolves continuously from one phase to the next. The 
complete mode function during all three epochs is 



00/, k) 



<f>£\v,k) 



for —00 < 77 < r]i deSitter, 



ocTad^itdiV, k) + Aad^ack 7 ?, k ) for f]i < i] < i] 2 radiation, 

. « fimltiv, k)+fl 0mat(^, k) for 7] 2 < 77 matter, 
where the coefficients a and (3 are given by 



(4.17) 



a (3 
(3* a* 



a (3 
13* a* 



rad 



a (3 
f3* a* 



(4.18) 



The mode function ( 4.171) is the normalized, continuous, graviton mode function which 
appears in the expression for the correlation function ( J2.65 ) . This expression for the mode 
function is exact, and valid for all wavenumbers k. 
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B. Corrections to the Instantaneous Phase 
Transition Approximation 



Our inflationary cosmological model undergoes instantaneous phase transitions; first be- 
tween the de Sitter and radiation phase, and then between the radiation and matter phase. 
At these transitions, the scalar curvature of the universe changes abruptly, since the second 
derivative of the scale factor ( |4. 1| ) is discontinuous. This abrupt change in curvature pro- 
duces gravitons, much in the same way as an abrupt change in the electromagnetic potential 
produces photons. This instantaneous phase transition is a good approximation, except at 
high frequencies, where it predicts too much graviton production P9fl . 

The physical universe transforms smoothly from phase to phase, with each transition 
taking place during a characteristic period of time. If the characteristic time of a phase 
transition is At, then one would expect the spectrum of gravitons produced by the phase 
transition to be supressed above a cut-off frequency f cnt , with / cut ~ 1/At. Equivalently, 
the production of gravitons whose wavelength is less than A cut = At is supressed. The ideal- 
ization that the phase transitions are instantaneous is a good approximation for frequencies 
below / cut . For this reason, the multipole moments (af) for small values of / should be 
unaffected by this idealization. The (af) for large I, however, will be overestimated if we do 
not "smooth out" the phase transitions. 

The adiabatic theorem [|TJ|,f)(| provides a simple way to account for the effects of "smooth- 
ing out" the phase transition, which does not require any detailed information about how 
the abrupt change in d(rj) is smoothed. The cut-off wavelength A cut corresponds to a cut-off 
wavenumber k cut - The adiabatic theorem implies that the Bogolubov coefficient (3 in ( [4.1 7| ) 



whose modulus squared gives the number of gravitons produced in a given mode |2J| , should 



decay exponentially for k > fc cut , while the Bogolubov coefficient a goes exponentially to I. 
So for large values of k with k ^> k cut no graviton production takes place. 

To be specific, consider the transition from the radiation to the matter phase. We will 
assume that the characteristic time of the phase transition equals the Hubble length at that 
time. This gives a cut-off wavelength 

AcutM = ^ = H-\V2). (4.19) 

At later times, the cut-off wavelength is redshifted by the cosmological expansion to the 
longer wavelength 

A C ut(r/) = A cut (r ?2 )44- ( 4 - 20 ) 
a{r] 2 ) 

The cut-off (comoving) wavenumber is then given by 

, cul = = a^-i. (4.2i) 

Acut W 

The adiabatic theorem now implies that for comoving wavenumbers above k cnt 



(k> k 

cut ) 



l + |/Wfc cu t)| 2 e 2 - 2fc / fc - 



1/2 



CWt^cut), (4.22) 



1 + |/Wfccut)| 2 

/W* > fc cut) = Anat^cuOe 1 -^. (4-23) 
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We use these formulae to determine a mat and /3 mat for k > k cut . They only significantly effect 
multipole moments (af) with I •> 1000. A similar analysis shows that the instantaneous 
transition from the de Sitter to the radiation phase only affects the moments with extremely 
large /. 



C. The Long Wavelength Approximation 

As noted in the introduction, the previously published calculations determine the angular 
correlation function using a "long wavelength" approximation to the graviton mode function. 
(We assume that the last-scattering event took place after the universe became matter- 
dominated, ie. r]2 < rj e ; for the rest of this paper, the "mode function" means the mode 
function during the matter phase). The long wavelength approximation is the same as an 
approximation for small wavenumber k. To make a small k approximation to the mode 
function, it is helpful to express the mode function in terms of spherical Bessel functions. 
Using (|4.7p , Q4.10D , and ( |4.17| ), one can write the mode function as 



0(7?, k) = -4i J±£* *!3L!»J [a + fljx{kin + r) 2 )) -i[a- (3] Vl (k( V + V2 )) \. (4.24) 



3tt p p (77 + 7/2) 

To understand the small k behavior of (f>(rj, k), one can expand the combinations of Bogol- 
ubov coefficients in square brackets as power series in k. One finds 

|a + 31 = T^ + °< r!) 



V 45?7i 

Furthermore, the small k behavior of the spherical Bessel functions can be understood by 
noting that for small argument, 

z < 1 =$> ji(z) ~ - and y x (z) ~ -z' 2 . (4.26) 
3 



Using ( |4.25| ) and (|4.26|) , and noting that 



3i(z) = y|;-W(*), (4.27) 
one sees that for small wavenumber k one may approximate the graviton mode function by 

0(77, k) « ./r2^ for V2 < V (matter). (4.28) 

The validity of this small k or long wavelength approximation, in the context of the angular 
correlation function, is discussed in section |VT| . 

Using ( 4.24 ) and (|4.25|) one can see why the long wavelength tensor perturbations can 



be thought of as classical gravitational waves. The Bogolubov coefficients are restricted by 
the constraint |a| 2 — |/5| 2 = 1, and in realistic inflationary models the "occupation number" 



\f3 1 [56] is much greater than one for small wavenumber k. Hence for small wavenumber a 



and (3 are both very large, and almost equal, as is apparent from ( [4.25 ). In this limit the 



graviton mode function (|4.28|) may be thought of and treated as a classical gravitational 



wave, as one would expect for a bosonic field with large occupation number. 
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D. Standard Results 



Using the long wavelength approximation to the graviton mode function we can repro- 
duce the standard results [f2~3|-|2~7|| for the angular correlation function due to gravitational- 
wave perturbations. Although the long wavelength approximation is valid only for small 
wavenumber k, we assume it to hold for all k. Recall that the first derivative of the mode 
function ( 2.40 ) appears in the angular correlation function. With the approximate mode 
function (4.28) and the definition ( |2.4U| ) for the function -F(A, k), one finds 



F(X,k) 



/ 12 Pds J5/2(k(r] 2 + r]e + A)) 
p P A; 3 / 2 (r/ 2 + r/ e + A) 3 / 2 " 



(4.29) 



Note that we have used the standard recurrence relations (equation (9.1.27) of reference 
|57| ) for Bessel functions to put -F(A, k) is this form. Substituting this into (|2.65|) and using 
( |4.27| ) one obtains 



6vr^ k~ 3 



dX 



J5/2(k(V2 + Vc + A)) Jl+l/2(k(llobs -Ve- A)) 



p P J Q (r)2 + Ve + X) 3 / 2 (r] ohs - rj e - A) 5 / 2 

To compare our results with standard formulae, define dimensionless variables 

x = k(i] 2 + i] c + A), 
b = k(i] ohs + r) 2 ). 

In terms of these variables, the multipole moments ( |2.67| ) are given by 

<k_r2 

y 



(4.30) 



(a 2 ) = 24vr 3 ^(/-l)/(/ + l)(/ + 2) f 
P P Jo 



where 



dx- 



Jb/2(x) Jl 



1+1/2 



(y 



X 



-'V 



X 



3/2 



(y — x) b / 2 



and 



V2 + Ve 

Vobs + m 



:i + z h 



'1/2 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



Equations ( |4.32j ) and ( |4.33| ) are equivalent to equation (8) in reference [26]. The lower limit 
of the integral in ( [4.33| ) appears different since the conformal time in our scale factor ( (4.1| ) 
during the matter phase is shifted from that in [26| by the constant rj 2 . 



V. PROGRESS TOWARDS A CLOSED FORM 
FOR THE ANGULAR CORRELATION FUNCTION C( 7 ) 

Equation ( |1.3| ) is a closed form for the angular correlation function due to scalar per- 
turbations. In this section, we attempt to find a closed form for the angular correlation 
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function due to gravitational- wave perturbations. Using (|2.65|) , (|2767|) , and ( |2.69|) one may 
write the correlation function as 







where 

B(r,^A;,7)=E( 2i + 1 )773^^( fcr )^(^ , )^(cos7), (5-2) 

and -D(A) is defined in ( 2.11|) . To obtain a closed form for the correlation function one must 
complete the integrals over A, A', and k, and the infinite sum over /. 



A. The Sum Over I 

One can sum over I and find a closed form for B(r,r' , fc, 7) using an addition theorem 
for spherical Bessel functions. Consider the addition theorem (see equation 10.1.45 in |57| ) 



^ = £(2/ + 1)^)^)^087), (5-3) 
where the length s is defined by the non-negative root of 

2 2 /2 / / \ 

s =r +r — 2rr cos 7. (5-4) 

The right-hand side of ( |5.3| ) is the same as the right-hand side of ( p.2| ), apart from the ratio 

of factorials. The ratio of factorials is just a fourth order polynomial in I. To generate this 
polynomial we define the derivative operator 

^ 1 9 . d 

V=- — — shi7— , (5.5) 
sm 7 (77 c/7 

which is the Laplacian on the unit two-sphere for functions with azimuthal symmetry. The 
Legendre polynomials are eigenfunctions of V, and obey 

VPi(cos-f) = -1(1 + l)Pi(cos7). (5.6) 

Using this one can quickly show that 

V(V + 2)P ; (cos 7 ) = |^-||pKcos7). (5.7) 

Using the addition theorem ( |5.3| ) and ( |5.7| ) one obtains 

B{r, r', k, 7) = V{V + 2)^^. (5.8) 

ks 
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One may distribute the derivative operator V(V + 2) on sin ks/ks to obtain a closed form 
for B(r,r',k,j). We prefer not to distribute the derivative operator, and instead use ( |5.1|) 
and (15.80 to write the correlation function as 



° h) = 7T I dX I dX L Tk*D*(\)D*w v{v + 2) ks(x,x^) ' (5 - 9) 

where now 

s(A, A', 7) = ^Jd 2 {X) + D 2 {X') - 2D{X)D{X>) cos 7 . (5.10) 

This form of the correlation function is very general. It only depends upon the cosmological 
model through the graviton mode function (or more precisely, its first derivative) which 
appears as F(X, k). 



B. The Integral Over k 

The next step to finding a closed form expression for the correlation function is to evaluate 
the integral over the wavenumber k. Since the derivative of the graviton mode function 
depends on the wavenumber k, one can not integrate over k without using a specific form 
for F(X,k). We use the long wavelength approximation ( |4.29|) , and assume it valid for all 



wavenumbers k. The accuracy of this assumption is discussed in section [VT]; the principle 
conclusion is that C{j) will be accurate for 7 greater than a few degrees. Substituting the 
long wavelength approximation into ( |5.9|) one obtains 



lobs"* Ijobs-* 

CM = 12^ [ dX [ dX'{D 2 (X)D 2 (X / )R(X)R(X')}- 1 
p p J J 





00 



where 



dk j 2 (kR(X))j 2 (kR(X')) | BmkajXA'^) (5 n) 

Jo k k A ks(X,X', r y) 



R(X) = 772 + Ve + A. (5.12) 



Note that the only dependence of the right-hand side on 7 is in the derivative operator 
V(V + 2), and in s(A,A',7). 

The derivative operator V{V + 2) is independent of k, so one might wish to take it outside 
the fe-integral. The remaining integrand could then be recast as a sum of trigonometric 
functions times powers of k. The problem with this is that the resulting integral over k is 
logarithmically divergent because the remaining integrand diverges as k~ l for small k. 

Still, one may take the derivative operator outside the integral by setting the lower limit 
to a small, positive number e. After applying the operator V(V + 2) one can then take the 
limit as e vanishes. So one can write the correlation function as 



fob S -1c Vobs-Vc 



C( 7 ) = 12^ J dX J dX'{D 2 (X)D 2 (X')R(X)R(X')}- 1 
P p 

xlim V(V + 2)K e (R(X),R(X , ),s(X,X',j)), (5.13) 
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where 



KfabX-r)) = J e dk _. (5.14) 

The function K e is well defined and finite for e > 0; one may evaluate it using standard 
techniques. 

To evaluate K e , express the spherical Bessel functions as exponential functions divided 
by powers [[54], expand the integrand, and integrate term by term (see 2.324.2 of reference 
|58fl). This yields 

K e {a, b, c( 7 )) = - a — In e + U{a, b) + V(a, b, c( 7 )) + 0(e), (5.15) 

where the functions U and V are independent of e, and terms which vanish as e goes to zero 
are not explicitly shown. The term proportional to lne and U(a, b) do not depend on 7 and 
are annhilated by V(V + 2). Only the function V(a, b, c( 7 )) contributes to the correlation 
function C( 7 ). The function V(a, 6, c( 7 )) is a sum of more than 25 terms, each of which is a 
rational function of a, b, and c, or a rational function of a, b, and c times In \p(a, b,c)\, where 
p(a, b, c) is a second order polynomial in a, b, and c. Using (|5.13|) and ( 5.15|) one can write 
the angular correlation function as 

c(7) = 12 * / / rfAP(p + 2) ^llgj ' (5 ' 16) 

where we have neglected to explicitly write out the function V, since it is not very illumi- 
nating. 



C. The Integrals Over A and A' 

The final step in finding a closed form for the angular correlation function is to complete 
the remaining integrals over A and A'. These integrals, however, are difficult for a number 
of reasons. Distributing the derivative operator V(V + 2) over the integrand of (|5.16|) yields 



on the order of 1000 terms. A large number of these terms are proportional to a logarithm, 
with argument linear in the function s(A, A'). The function s(A, A') is the square root of 
a second-order polynomial in A and A' which is not factorable for arbitrary 7 . Integrating 
terms like these over A and A' is not trivial. Other terms are proportional to odd powers of 
s(A, A'), and are difficult for the same reason. The total number of terms, combined with 
the difficulty of integrating each term, impedes further progress. 

Other methods for finding a closed form for the angular correlation function do not appear 
more promising. One can write K e in the limit as e vanishes as a hypergeometric function 



of two variables (see 6.578.1 in reference ||58||), but again the remaining integrals over A and 
A' are difficult. The integral over the wavenumber k can be evaluated before summing over 
I, though this involves the integral of four Bessel functions, each with a different argument. 
One may also consider the integrals over A and A' first. These integrals are almost, but 
not quite, standard integral transforms of Bessel functions. Another approach is to begin 



with the Sachs- Wolfe operator ( 2.37| ) and calculate the angular correlation function ( J2.42 ) 



directly without any expansions in terms of spherical harmonics. This approach, however, 
reproduces (|5.9| ). 
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VI. COMPARISON OF THE EXACT AND LONG WAVELENGTH 
APPROXIMATE MULTIPOLE MOMENTS (af) 



A. Analytical Comparison 

We are considering an inflationary cosmological model that begins with a de Sitter phase 
followed by radiation- and then matter-dominated phases. The graviton mode function 
( |4. 17] ) that we obtained for this model is exact, and valid for all wavenumbers k. Using this 



mode function one can calculate the multipole moments (af) an observer in this universe 
model would measure. 

In the standard literature p3|- [26|j2^j38|j4^ -pi| , however, the long wavelength approximate 



mode function ( [4.28 ), rather than the exact mode function ( 4.17 ), is used to calculate the 
multipole moments. The approximate mode function is only valid for longer wavelengths, 
so we expect the angular correlation function, when calculated using the approximate mode 
function, to be accurate only on large angular scales. Equivalently, we expect the moments 
(af) calculated using the long wavelength approximate mode function to only be accurate 
for small I. 

The straightforward way to determine for which (af) the long wavelength approximation 
is valid is to numerically calculate the moments using both the exact and approximate mode 
functions, and compare. To calculate the moments one substitutes into ( |2.65| ) either the 



derivative ( [4.29| ) of the approximate mode function, or the derivative of the exact mode 



function fl4.17| ). One then uses Q2.671 ) to obtain the multipole moments. Making these 
substitutions, one finds for the moments calculated with the long wavelength approximate 
mode function 

2 2 (I + 2)! pds f°° 3 ~2 

(fll ) long wavelength approximation 487T — — r— / dy y (|/), (6.1) 

[l- 2)1 p p Jo 

and for the moments calculated with the exact mode function 

„2\ ylQ 2 + 2)' Pds f°° , 3 fry, / \x2/„.\ , "T f . .\\r2 I 



(/ - 2)! p p Jo 

The dimensionless variables x and y are defined by the change of variables 

x = — — — — - and y = k(r] ohs - rj e ). (6.3) 

Vobs - Ve 

The functions Ji(y) and Yi(y) are defined by 

j,(y) = f 1 dx J wMyx)j 2 (y^-x)) ~ [i Ji +1/2 {yx) y 2 m-x)) 

where £ is a dimensionless constant determined by the redshift of the last scattering surface 

£ = r]ohs + 7]2 = [1 - (1 + Z Xs )- 1 ' 2 }- 1 . (6.5) 
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In realistic cosmological models, £ is slightly greater than one. The three functions Tj(y) 
depend on the Bogolubov coefficients a and (3, and are defined as 

2 



Ti(y) = 
T 2 (y) = 



y d CfCa(a + /9) 



2/ 3 Ci 2 C 2 (a-^) 



T 3 (y) = ^ 6 a 2 Im{a^}. 



(6.6) 
(6.7) 
(6.8) 



The Bogolubov coefficients are given by (|4.18|) with the changes of variable (|6.3| ). The 
dimensionless constants Ci an d C2 are determined by the redshifts Z end , Z equa i, and Z\ s 
defined in (gg-gg): 



Ci 

C2 



1/ yr+zi 



l + Z 



equal 



r] obs - r] e 2 V V 1 + z \s - 1/ (1 + ^cnd) 

1 



'/2 



1/ v / T+^ 



r/ obs -r/ c 2V v / l + ^is-l/ + 



(6.9) 
(6.10) 



J equal 



To determine for which (af) the long wavelength approximation is valid one numerically 
integrates ( |6.1| ) and ( |6.2| ) and compares the values for the moments. 

For cosmological models with "enough" inflation to solve the horizon and flatness prob- 
lems, £1 is very small since Z end > 10 26 . For this reason one can approximate the Tj(y) 
by 

(43/C2 cos(yC 2 ) - sin(yC 2 ) + 8y 2 (l sin(yC 2 ) + sin(3yC 2 )) 2 



Tx(y) 

T 2 (y) 
T 3 (y) 



36y 2 C 2 2 



+ 0(Ci), 



%C2 sin(yC 2 ) - cos(yC2) + 8y 2 C| cos(yC2) + cos(3yC 2 ))' 

362/ 2 Cl 



+ 0(d 



(6.11) 
(6.12) 



VCr 2 



sin 2 ( 2 /C 2 ) sin(%C 2 ) - 4yC 2 (l + 2 cos(2yC 2 )) sin%C 2 ) - 32y^2 cos(yC 2 



+16y d C| cos(2j/C2) + 16y 4 C 2 4 sin(2yC 2 ) + O(Ci). 



(6.13) 



In what follows, we neglect the 0(d) and higher terms in Tj(y). Note that the standard long 
wavelength approximation (|6.1|) is equivalent to setting Ti(y) = 1 and T 2 (y) = T 3 (y) = 
in the exact expression (|6.2|). Indeed, expanding ( p. lip , ( |6.12j ), and ( |6.13| ) as power series 
in y one finds 

T 1 ( 2/ ) = l + 0(yC2) 2 , 



T a (y) 



256 



18225 



(yC 2 ) w + 0(y(2) 



12 



32 
135 



(2/C2) 5 + 0(y( 2 ) 7 . 



(6.14) 
(6.15) 
(6.16) 
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So for y( 2 < 1, to a good approximation one has Tx(y) = 1 and T 2 (y) = T 3 (?/) = 0. 

Using the power series Q6.14 - B.16 ) we can understand why the standard approximation 
( |6.1| ) is the long wavelength approximation to (p.2|). The functions Ji(y) and Yi(y) are 
peaked near y = I. Figure |l| shows Ji(y) and Yi(y) for I = 10 and I = 100. Hence, if I < Q 2 , 
Ji(y) and Yi(y) only have support for y < ( 2 ■> which is the same range for which Ti(y) ~ 1 
and T 2 (y) ~ T 3 (?/) « 0. For y > Q 2 > an d have no support, and the second and 
third terms in the integrand of the exact formula fl6.2| ) do not contribute for large y. So 
one expects the standard approximation ( |6.1| ) to give accurate values of (af) for I < (^T . 
From ( |6.10|) note that for realistic models, ( 2 ~ ^equal- ^° one ex P ec ts that for I < the 
standard long wavelength approximation gives accurate values for the multipole moments 
(af)- 



B. Numerical Comparison 

Table f| lists the multipole moments for various I values, calculated using both the 
approximate formula (|6.1|) and the exact formula (|6.2|) for Z en d = 10 27 , Z equa i = 10 4 , and 



Z]s = 1300. The difference between the exact and long wavelength approximate moments is 
shown in Figures ^ [| for different values of the cosmological parameters. For 2 < / < 10 our 
values of (af) agree very well with those of White |^3| (due to a difference in the definition 
of (af), our results are smaller than White's by a factor of 2/ + 1). As expected, for smaller I 
the values of (af) from the approximate formula are in good agreement with the exact (af). 
For I < 30, the difference between the exact and approximate moments is less than two 
percent of the exact result. For I < 100, the difference is less than twenty percent. When / is 
200, however, the disagreement is more substantial; the exact value is more than twice the 
approximate value. The disagreement is even more for larger /, and for I = 1000, the exact 
value is a factor of 69 larger than the approximate value. The long wavelength approximate 
formula (|6.1|) substantially underestimates the contribution of the large / moments (af) to 
the angular correlation function C ('-/). 

Similar results, which reveal the shortcomings of the approximate formula for (af) have 
been obtained by Turner, White and Lidsey |53 |. Their approach is less analytical than 



our own; they use numerical methods to solve the Klein-Gordon equation and obtain exact 
mode functions 0(?7, k) analogous to our equation ( |4.17|) . They express these exact solutions 



in terms of the standard long-wavelength approximate mode functions, using a "transfer 
function" . Figure [| shows the results of our best attempt to obtain the Turner, White, and 
Lidsey results from our analytical formula, together with their published data. By tuning the 
parameters of our cosmological model to Z\ s = 900 and Z equal = 2500, we have been able to 
obtain fairly close agreement between the two sets of results. One should note, however, that 
Turner, White, and Lidsey consider a universe which is not completely matter-dominated 
at the time when the CBR is emitted. Their universe model is more realistic than our own, 
since we have assumed the universe to be completely matter-dominated at the time of last 
scattering. Ng and Speliotopoulos |[59|| , using a WKB formalism, have also considered a 
universe which transforms smoothly from the radiation- to the matter-dominated phases. 



Some of their results, however, do not appear consistent with those of |53 . 

We have shown that the long wavelength approximate formula ( |6.1| ) substantially under- 
estimates the contribution of the large / multipole moments (af) to the angular correlation 
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function C{pf). Although this long wavelength approximation has been used previously to 
interpret published experimental data, one does not expect the new results presented here 
to significantly affect the conclusions. This is because for reasonable values of the redshift 
Z e qaah the discrepency between the approximate and exact results is significant only for 
multipole moments which one expects would be dominated by the contribution from scalar 
perturbations 



45| . However, Krauss and White |42] and Grishchuk have suggested 



that the relative contribution to the CBR anisotropy from gravitational waves has been 
underestimated, and that these contributions might dominate the multipole moments. 



VII. CONCLUSION 

In this paper, we have shown how the rapid expansion of the universe during an inflation- 
ary phase creates large numbers of gravitons, whose collective effects produce potentially- 
observable fluctuations in the temperature of the CBR. The correlation function of these 
temperature fluctuations may be calculated from first principles; for example the overall 
magnitude of the perturbations is determined by the uncertainty principle. The exact ex- 
pression that we obtain for the correlation function agrees with standard published results 
for the lower multipole moments, but has larger temperature fluctuations in the higher mul- 
tipole moments than predicted by the standard published formulae. This appears to be in 
good quantitative agreement with recently published numerical work by Turner, White, and 
Lidsey |kJ . The larger predicted temperature fluctuations in the higher multipole moments, 
however, most likely will not lead to a reinterpretation of the experimentally observed data 
since it is generally expected that the observed anisotropy for the higher multipole moments 
will be due almost entirely to scalar, rather than tensor perturbations. 

As mentioned in the introduction, the original discovery that a rapidly-expanding uni- 
verse could create relic gravitational waves was made by Grishchuk p0| . In recent work 



p8| . |49|| , Grishchuk analyzed the temperature fluctuations produced by these waves, using 
the techniques of quantum optics. In his analysis, the classical gravitational field "inter- 
acts" with the gravitons and acts as a "pumping" field. This leaves the gravitational field in 
a squeezed quantum state today. Grishchuk stresses the importance of the resulting phase 
correlations to the final form of C(j). 

In our language, the (quantized) gravitational field is taken to be in the vacuum state 
of the initial de Sitter phase. (Note that we use the "Heisenberg picture" of quantum fields 
in which the states do not evolve with time, but the operators do; we also assume that a 
sucessful inflationary stage leaves the universe indistinguishably close to the de Sitter vacuum 
state). Although we do not use any of the techniques of non-linear quantum optics that 
Grishchuk advocates, we nevertheless reproduce, as intermediate results, his final formulae 
for C ('-/). In particular, Grishchuk's formula (10) from reference |0| is the same as our 
equation (|5.9|) with 7^0. Formula (11) from reference [HS[ is the same as our equations 



( p.43| , pT46| - p.48| ), and formulae (12-13) from reference (|8| is the same as our equation ( |5.9| 



with the action of the operator V(V + 2) expanded out. We agree that the correlation 
between phases is important; in the sense that for example in our equation ( |6.8|) the value 
of T 3 depends upon the relative phase of the positive and negative frequency wavefunctions. 
However, we stress that results identical to Grishchuk's may be obtained, as we have shown, 
using only the standard machinery of linearized quantum fields in curved spacetime [r5,5o|. 
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APPENDIX A: 



In section |IIE2| we argued that based on the isotropy of the initial state of the universe 
(which we took to be the de Sitter vacuum state), and on the istropy of the FRW model, 
one expects the angular correlation function to be rotationally invariant. For this reason 
one may write the matrix element {0\C^ q Ci m \0) as in fl2.44| ), and then use fl2.41| ) for Ci m to 
solve for (af). In this appendix we sketch this calculation. A somewhat more complicated 



version of this calculation may be found in |43 



The primary advantage of writing the matrix element as in ( |2.44j ) is that it allows one to 
make a useful choice of coordinates and evaluate the integrals over angular variables. Using 



(p.41|) for Ci m one obtains from (|2.44|) 



<«?> = lJ dX 'J dX J —F(X',k')F*(X,k) 



e ab (k e K d (k e ) + e* ab Ck e )e cd (k e 



A c A c 

x f dfl% f ^^(v^^^v^^^e-^^W^-- ^'^), 



(Al) 



where we have set p — I and q = m to eliminate the Kronecker delta functions on the right- 
hand side of ( |2.44|) . Since by assumption both sides of the equation above are independent 
of m, one may sum both sides from m = — I to m = I. Using the addition theorem for 



spherical harmonics [53], and cancelling factors of (21 + 1) on both sides, one obtains 

-^obs 



(«?> 



16n 



x 



dX' J dX 

A c A c 

dVL 



d 3 k r - 

— F(X', k')F*(X, k) \e ab (k e )e* cd (k e ) + e* ab (k e )e cd (k £ 



where the angle 7 is defined by 



cos 7 



U V n 



(A2) 



(A3) 



Note that we have not yet made a specific choice of coordinates. 

One is free to choose whatever coordinates one wants to compute the integrals over the 
angular variables flu, fly, and fl^. In particular the choice of coordinates for fl u and fl v may 
depend on the vector k c . We choose coordinates so that the vectors u c and v c are written 
in terms of the (rff, n c , k c ) triad as 



u c = sin 8 U cos 4> u m a + sin 9 U sin 4> u n a + cos 9 u k a 
v c = sin 9y cos <fi,~ u m a + sin 9% sin n a + cos 9 v k . 



(A4) 
(A5) 
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With this choice of coordinates, and using the form of the polarization tensors given in 
C 2.26| ) , one can quickly show that the contraction between the polarization tensors and the 
unit vectors is 



e ab (k e )e: d Ck e ) + e: b (k e )e cd (k e 



v°v d u a u = sin 2 6u sin 2 9$ cos(20^ 



(A6) 



Also using these coordinates one may again use the addition theorem for spherical harmonics 
to write 



47T - 

P;(C0S7) = Zl Y lm( u,^u)Ylm(Ov,(t>v), 



(A7) 



where 



cos 7 = cos Ou cos 9~ + sin 9% sin 9% cos(0£ 



(A8) 



(Note the following subtle point. The Y l ^ n (9{ i ,,fa) and Yi m (9u,fa) in the right-hand side 
of ( |A7p do not in general have the same values as the spherical harmonic functions which 
appear in ( |A1| ) because we have done a coordinate rotation that depends on k . In general 
these values are related by a linear expansion involving Clebsch- Gordon coefficients.) Using 
and (|A7|) one finds for the multipole moment 



-^obs -^obs 



( a f) = _L J d \> J dX J ^A F (X',k')F*(X,k) J da b J ^. e -^(A)cos^^(V)cos e ,) 

x sin 2 6 & sin 2 9, cos(20 fl - 2 fa) I ]T YfjO^ fa)Y lm {9,, fa) \. (A9) 

V ) m=—l 

In this form one may evaluate all the integrals over angular variables. 

The integrals over the angles 9u,9y,4>u, and fa can be done in a straightforward way by 
writing the spherical harmonics as products of exponentials and Legendre functions. With 
the definition 



Y, 



(21 + I) (I - m)\ 
\| 4vr (l + m)\ 



P] n {cos9)e 



im<f> 



(A10) 



one finds after integrating by parts I — 2 times (formula 3.387.2 in Gradshteyn and Ryzhik 
|H is helpful) 



2 _ zl ii±n 

{ai) ~ 2 (Z - 2)! 

o k 



X 



^obs 


Aobs 






A e 


A e 


sin % j 




' % 



7i(fcD(A))j,(fcD(A')) 
' k 4 D 2 (\)D 2 {\>) 



^ (^m,2 + ^m,-2)- 



(All) 



m=— I 



The remaining integals over 9^ and are trivial and yield 4tt. The sum over m is also 
trivial and contributes a factor of 2. So one obtains for the multipole moment 



35 



A c A c 

Recalling the definitions of D(X) and I\{k), one can write this as 



This result is the same given in (|2.67|) . 



APPENDIX B: 

This appendix describes the numerical techniques used in section [V|. 

The primary numerical technique used to evaluate both the approximate ( |6.1| ) and exact 
multipole moments ( |6.2| ) is numerical integration. Both integrals over y in (|6.1|) and ( |6.2|) 
were done using a fifth order embedded Runge-Kutta-Fehlberg algorithm with adaptive 



stepsize control [61]. Although formally the upper limit of the integral extends to infinity, 
we only integrated until the remaining contribution became negligible. This is possible 
because the integrands in (|6.1|) and (|6.2|) fall off at least as fast as y~ 2 for large y. Special 
care must be taken in determining when the remaining contribution is negligible since the 
integrand does have periodic zeroes, even for large y. 

Both Ji(y) and Yi(y) (or more precisely, these functions multiplied by y 7 ^ 2 ) were also 
calculated using a fifth order embedded Runge-Kutta-Fehlberg algorithm with adaptive 
stepsize control. No special treatment is needed since both the upper and lower limits are 
finite, and the integrands are well behaved. The spherical Bessel functions in the integrands 
of ( |6.4|) can be expressed in terms of trig functions |54|] , and evaluated using standard 
machine routines. 

The Bessel function with index 1 + 1/2 was evaluated with the routine "bessjy" given 
in chapter 6 of Although this routine is very accurate and fairly fast, we did not use 



it to calculate the value of the Bessel function every time it was needed in the integration 
algorithms. This is because the argument of the Bessel function is a function of both x 
and y, and so the argument is unique for every step taken while the integral over x is 
being calculated; it is not possible to store certain values and "reuse" them later. A typical 
integration to find a single moment for a particular I would easily require on the order of 
10 6 calls to the routine bessjy. 

To reduce the number of "expensive" calls to bessjy we used a cubic spline interpolation 
scheme to calculate the Bessel functions with index / + 1/2. For each different value of 
a table of Bessel functions evaluated at equally spaced intervals A = 1/32 was tabulated 
using bessjy. This interval is small enough so that cubic spline interpolation gives values 
accurate to at least one part in 10 6 . The cubic spline was done using the routines "spline" 
and "splint" from chapter 3 of |H[ . 
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FIGURES 



FIG. 1. The functions \Ji(y)\ (solid curve) and (dashed curve) for I = 10 and 100, 

normalized so that their maximum value is one. For all curves Z\ s = 1300. Both Ji{y) and Yi(y) 
peak fairly strongly near y = I, which makes the long wavelength approximate formula for the 

1 /2 

angular correlation funcion multipole moments accurate for I < Z' ,. 



FIG. 2. Multipole moments (af) normalized to the quadrupole moment (a 2 ), with 
Mi = ( — ) ( ^t" 1 ) ] [ jlX j . The upper curve is the correct result calculated using the exact 



Pds J \ 6 J \(a%) 

graviton mode function. The lower curve is the result obtained from the standard formula found in 
the literature using a long wavelength approximation to the graviton mode function. Both curves 
have Z cnd > 10 20 , Z cqua i = 10 4 , and Z h = 1300. 



FIG. 3. Multipole moments (af) normalized to the quadrupole moment (a 2 .) with Mi the same 
as in Figure ||. All three curves are calculated using the exact graviton mode function, and have 
Z en d > 10 20 and Z equa i = 10 4 . The upper curve has Z\ s = 1300, the middle curve Z\ s = 800, and 
the lower curve Z\ s = 400. 



FIG. 4. Ratio of multipole moments obtained with the long wavelength approximation to the 
exact multipole moments with R/ = (a?), , /(a 2 )- All three curves have the same 

^ 1 \ t ' long wavelength approx / \ t / 

redshifts as in Figure ||. The approximate moments fail to be accurate for / > Z^\^. 



FIG. 5. Multipole coefficients (af) normalized to the quadrupole moment (a 2 .), with Mi the 
same as in Figure ||. The diamonds show the results of Turner, White, and Lidsey |53| obtained by 
expressing exact mode functions (obtained by numerically integrating the massless Klein-Gordon 
equation) in terms of the standard long-wavelength approximate mode functions using a "transfer 
function". The upper curve shows an exact result obtained from our analytic formula ( |6.2[ ) with 
Z ls = 900 and 

■^equai — 2500. These parameters were chosen because they appeared to give the 
best match to the Turner, White, and Lidsey result. The lower curve is the result obtained from 
the standard long wavelength formula ([O]) with the same parameters used for the upper curve. 
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TABLES 



TABLE I. Change of dependent and independent variable needed to cast the massless 
Klein-Gordon equation in the form of Bessel's differential equation, and positive frequency (unnor- 
malized) solution. 

Epoch Dependent Independent Solution <j) 



-oo < 77 < 7/1 cf> = (7/ - 2r] 1 ) 2 x 

771 < 7] < 772 4> = x 

772 < 7? = (7/ + 772)~ 1 X 



TABLE II. Multipole coefficients (of) for various / predicted for a stochastic background of 
gravitational radiation generated by exponential inflation. Exact values are calculated using the 
exact graviton mode function in ( |2.67| ) for the multipole moments. Approximate values are cal- 
culated using the standard long wavelength approximation to the graviton mode function. The 
values in this table are for redshifts Z en d > 10 20 , Z eq = 10 4 , and Z\ s = 1300. 

I Exact Approximate 



2 


1.55 


1.55 


3 


6.07X10- 1 


6.07X10- 1 


4 


3.44xl0 _1 


3.44xl0 _1 


5 


2.27xl0 _1 


2. 27x10^ 


6 


1. 62x10^ 


1. 62x10^ 


7 


1.23X10- 1 


1. 22x10^ 


8 


9.61xl0- 2 


9.59xl0" 2 


9 


7.75xl0~ 2 


7.73xl0~ 2 


10 


6.38xl0~ 2 


6.36xl0~ 2 


25 


1.08xl0~ 2 


1.07xl0- 2 


50 


2.11xl0~ 3 


2.01xl0~ 3 


75 


5.39xl0~ 4 


4.85xl0~ 4 


100 


1.22xl0~ 4 


LOlxlO" 4 


250 


6.45xl0~ 7 


1.47xl0~ 7 


500 


5.98xl0~ 8 


3.49xl0~ 9 


750 


1.49xl0~ 8 


4.93xl0~ 10 


1000 


5.07xl0~ 9 


7.33X10" 11 



z = k(r/ — 2r]i) 

z = kr) 
z = k{rj + 772) 



( V -2 m fh^(k( V -2 Vl )) 
(V + rj 2 y lh< i\k(r] + 772)) 
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